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Mimicking the description of spinning particles in General Relativity, the Fermat Principle is 
extended to spinning photons. Linearization of the resulting Papapetrou-Souriau type equations 
yields the semiclassical model used recently to derive the "Optical Hall Effect" for polarized light 
(alias the "Optical Magnus Effect"). 
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Light is an electromagnetic wave, whose propagation 
is described by Maxwell's theory. It can also be viewed, 
however, as a particle (a "photon"). Here we adhere to 
the second approach: we describe light by a bona fide 
mechanical model in that we use a Lagrangian. 

In traditional geometrical optics the spin degree of free- 
dom is neglected, and the light rays obey the Fermat 
Principle 0]. In the intermediate model advocated by 
Landau and Lifchitz 0, the photon is polarized, but the 
polarization is simply carried along by the light rays, and 
has no influence on the trajectory of light. Recent ap- 
proaches 0, 3 go one step further : the feedback from 
the polarization deviates the trajectory from that given 
by the Fermat Principle. A dramatic consequence is that, 
for polarized light, the Snel(-Descartes) law of refraction 
requires correction : the plane of the refracted (or re- 
flected) ray is shifted perpendicularly to that of the inci- 
dent ray [Jj. This "Hall Effect for light" is a manifesta- 
tion of the Magnus-type interaction between the refrac- 
tive medium and the photon's polarization Q . It can be 
derived in a semiclassical framework, which also includes 
a Berry-type term 0, El • 

In this Rapid Communication, we argue that the devi- 
ation of polarized light from the trajectory predicted by 
ordinary geometrical optics is indeed analogous to the 
deviation of a spinning particle from geodesic motion in 
General Relativity. The resulting equations are reminis- 
cent of those of Papapetrou and Souriau |8| ■ 

In detail, the Fermat Principle of geometrical optics 
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says that light in an isotropic medium of refractive index 
n = n(r) propagates along curves that minimize the opti- 
cal length. Light rays are hence geodesies of the "optical" 
metric <7y = n 2 (r)<5,y of 3-space. To extend this theory 
to spin we consider the bundle of positively oriented or- 
thonormal frames over a 3-manifold endowed with a Rie- 
mannian metric g^. At each point, such a "Dreibein" 
is given by three orthogonal vectors U l ,V l ,W l of unit 
length that span unit volume. We stress that the [6- 
dimensional] orthonormal frame bundle we are using here 
is a mere artifact that allows us to define a variational 
formalism. Eliminating unphysical degrees of freedom 
will leave us with 4 independent physical variables. 

Introducing the covariant exterior derivative associ- 
ated with the Levi-Civita connection, DU k — dU k + 
T k jdx i U 3 , we posit the reparametrization- invariant ac- 
tion 

f dx i f DW i 

S = SWmat + S'spin = K / U. t —dT~S / Vi ^ dr, (1) 

where r is some parameter along the light ray. The pa- 
rameters s and k > are interpreted as the spin and the 
color, respectively. Upon first quantization, k becomes 
indeed, for a monochromatic wave, 2ttH/X, where A is the 
wavelength For the photon s = ±h, but we keep it 
arbitrary for future convenience. Equation is supple- 
mented with the constraints UiU 1 = ViV 1 = WiW 1 = 1, 
and UiV* = U i W i = V i W i = 0. 

The first term in is [k times] the usual optical 
length; the second, "Wess-Zumino-type" ^(| term, that 
arises naturally in the geometric framework of , corre- 
sponds to the Berry connection, and is indeed analogous 
to the torsion term considered by Polyakov [l2j] . 

The Euler-Lagrange equations are obtained as follows. 
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Variation of the first term in yields 
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where <5rt^ fc = 5U k + T k m 5x l U m is the covariant variation 
of the vector field U % . For the spin term, straightforward 
calculation yields 
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where R jkU = ffi m (0,TJg - d k T" 



is the Riemann 



tensor and S lJ = — s(V l W J — W l V J ) is the spin tensor. 
Then the variational principle 5S — allows us to infer 
the pair of equations 
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where ll is a Lagrange multiplier which enforces the or- 
thogonality of U l and 5rU z . Inserting DU k /dr into 
Eqn Q and redefining the parameter along the ray by 
dt = (/j,/K,)dT yields the Papapetrou-Souriau type |8| 
equations, 
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where — Rij — \Rg%j is the Einstein tensor and the 



J %3 — 2 

matrix R(S) k — R^S 13 represents the interaction of 
spin with the curvature, responsible for tidal forces. Due 
to the spin-curvature coupling, the direction of the ve- 
locity differs, in general, from that of the spin vector 

sU l - - x T,r-> k s l ,. 

For the optical metric, the Christoffel symbols are 
T k j = ±(din6 k + djn6 k - d k n8 i: j). (The optical met- 
ric is hence not flat unless the refraction index is con- 
stant.) Putting r = (x l ), u % = nU 1 ^ 1 = nV l ,w l = nW 1 , 
introducing the momentum, 



p = n 
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and denoting the derivative w.r.t. t by a "dot", our La- 
grangian in Q can also be presented as L = p f — sv-w. 
The equations of motion for r and p read, in this case, 
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These equations describe spinning light in an inhomoge- 
neous medium. Let us mention, for completeness, that 
the evolution of the spin vector, which follows from JJJ, 
is given by su = -nc[r — ~ V(i) x r] x u. 

If the medium is spherically symmetric, n — n(r), 
conserved angular momentum is readily derived using 
Noether's theorem. It reads 



r x p + su. 



(12) 



Let us now discuss some particular cases of our general 
theory 

(i) For s — we have p 2 = n 2 K 2 , a = 1 and Ap = p. 
Introducing the elementary arc length da = riKdt, we 
recognize the usual Fermat equations, ndvjda = p/re, 
d{p/K)/da = Vn 0. 

(ii) In a homogeneous medium, n = const . we get, 
for any color, k, and spin, s, the same equations: light 
propagates along straight lines parallel to p = uku. The 
model is invariant w.r.t. the Euclidean group SE(3) con- 
sisting of space translations and rotations. The asso- 
ciated conserved quantities are the linear momentum, 
p = uku, and the angular momentum, (|12f> . which is 
now J = rxp| sp/p. 

(iii) In a medium with slowly varying refractive index, 
terms involving second-order derivatives and quadratic 
expressions in V(l/n) can be neglected, e.g., Rij — 
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light is approximately tangent to the spin and the latter 
is approximately parallel transported, 
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In p-terms, p 2 
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In the case of spherical symmetry, the general angular 
momentum (I12f) reduces, up to the approximately con- 
served extra term (s 2 /«;)V(l/n) x u, to the expression 
used by Onoda et al. in |3J , namely to 
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Let first consider the free case, n = 1. The variable r 
used so far has been an arbitrary point of the light ray. 
Now, the ray itself can be labeled by its direction u and 
q = r — (u ■ r) u, which is in fact the shortest vector 
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drawn from the origin to the ray (orthogonal to the unit 
vector u) , and can be thought of as the "position" of the 
ray. The 4-dimensional manifold, A4, of light rays de- 
scribed by u and q, has the topology of the tangent bun- 
dle of the two-sphere and can be identified with a coad- 
joint orbit of SE(3). The Casimir invariants of the orbit 
(which determine unitary irreducible representations) are 
k = p (= vP 2 ") and s = 3 -p/p. The corresponding orbit, 
M, is endowed with the canonical symplectic structure 

uiQ = ndui A dq l — —e^k^du^ A du k , (16) 

see 9] . The monopole-like term in (II 61) is the Berry cur- 
vature, TjSCij^dpJ A dp k /p 3 . It makes the components of 
q non-commuting [JjJ, i- e -> Cartesian coordinates q\ and 
92 have non- vanishing Poisson bracket, {91,(72} = s/n 2 . 
Upon (first) quantization, in the case s = h, the quantum 
"position" operators 91 and 92 satisfy 

[91,92] =i(A/27r) 2 , (17) 

where A = 2ttH/k. The Heisenberg uncertainty relation 
read, therefore, A91 • A92 > ^(A/27r) 2 , which provide a 
new interpretation of the localization defect of spinning 
light rays, limiting the resolving power of optical instru- 
ments to the order of the wavelength, Aq sa O(X). 

In a non-trivial refractive medium the (exact) two-form 
(|T1)|> is replaced by 

w = nDUi A dx l - ^R(S) lj dx t A dx j 

-^V9^ 3 kU l DW ADU k (18) 

on the orthonormal frame bundle. The Euler-Lagrange 
equations ©-0 correspond in fact to the kernel of the 
two- form . see [ll| . Conversely, the spin term in our 
Lagrangian comes from a potential for the spin terms in 
the two-form (|18|1 . 

Now putting p — > p/n and t -4 «r, our linearized 
equations ifTlJ) become those proposed in £|. 

The relation to the model of Onoda et al. |3| is more 
subtle. In their approach, polarization is an additional 
variable, represented by a two-component complex vec- 
tor, z — (z a ) with a = ±, such that \z + \ 2 + \z^\ 2 = 1, 
acted upon by su(2). Their semiclassical equations of 
motion can be written as 

r' = — k + k' x n ab z a z b , k' = -V(-)k. (19) 
nk n 

supplemented with 

4 = fcV(i)-A ah z b , (20) 

where k is the wave vector and A afc (k) is an su(2)- 
valued non-abelian "Berry" vector potential. The Berry 
curvature can be represented by an su(2)-valued vector 
fi = <73k/fc 3 , a Dirac monopole in k-space, diagonally 



embedded into su(2). The vector potential for 17 can, 
therefore, be chosen as A ah = iA(a^) ab where A(k) is 
a monopole potential, rot A = k/fc 3 . The number of 
equations in l|19l) - (|20l) can be reduced to two. The po- 
larization equation l|20l) can in fact be solved formally by 
parallel transport, z a = e la6 z®, where the phase is given 
by the non-integrable phase factor 9 = Jk (V(^) • A)dcr. 
Then the Berry term becomes simply fl ab z a Zb = sk/fc 3 
where s — \z + \ 2 — |z-| 2 = |z+| 2 — \z°_ | 2 , since the \z a \ 2 
are separately conserved. Notice that this s is a constant 
of the motion, which can take any value between —1 and 
+1. Identifying the wave vector, k, with our momentum, 
p, and putting (•)' = (n 2 n)d/dt transforms finally H19|) 
into our equations (|14H . 

We did not consider the polarization in our framework. 
As long as we are only interested in describing light rays, 
polarization is a secondary quantity, whose only role is to 
generate spin, which in turn deviates the trajectory from 
that of conventional geometrical optics. It is hence more 
appropriate to speak of spinning light than of polarized 
light. Let us nevertheless mention that first quantization 
along the lines of Q of the classical model (|16|l . with 
Casimirs k and s = ±fi, yields, in the gauge divA = 0, 
the vectorial Helmholtz equation 

(A + fc 2 )A = 0, (21) 

where k = k/H. It follows that E = ikA and B = rot A 
satisfy the field equations 

rotE-ifcB = 0, (22) 
rotB + i/sE = 0, (23) 

associated with our Euclidean model. Promoting the 
parameter t as "time" , these are indeed the vacuum 
Maxwell equations for the stationary fields ~Ee~ tkt and 
Be _lfct , respectively (c = 1). In an isotropic medium, 
Eqn (|2*5|l is generalized Q to 

rotH + ifcD = 0, (24) 

where D = eE and B = /iH, with e and fi the permittiv- 
ity and permeability, respectively. The refractive index 
is n = y/ejl. Remarkably, the field equations can, again, 
be rewritten in terms of "optical" metric, namely in the 
form (I22|) and (|24|l above, replacing the operator rot 
by its curved-space form, rotX — * n~ 3 rot(n 2 X), and 
rescaling the fields, E — > n~ 2 E, H — > n~ 2 H, B — > ri~ 3 B, 
D — > n~ 3 D and e — > n~ 1 e, fi — » nT 1 [i. 

Conventional geometric optics can be derived from the 
eikonal approximation of Maxwell's electrodynamics 0] . 
Here we followed the opposite way : we started with a 
classical model and derived the stationary Maxwell equa- 
tions H22I24|I by (first) quantization. Although we have 
not yet been able to deduce our action J^l from taking 
a suitable semiclassical limit of I|22l24ll . we emphasise 
that our model actually comes from first principles - but 
of those of Mechanics. ^lj. Firstly, the free model is 
constructed along the lines suggested by Souriau's 0, 



4 



applied to the Euclidean group. The second step is min- 
imal gravitational coupling, which amounts to replacing 
the ordinary scalar product by the one associated with 
the "optical" metric and the ordinary derivatives by co- 
variant derivatives; this yields the two- form (|18fl . The 
latter is in turn associated with first-order variational 
calculus on "phase space" , whose Lagrangian is precisely 
our (JTJ). 

Our model reproduces, at first order, the phenomeno- 
logical descriptions proposed in 0,3 which can, in turn, 
be derived by taking an improved semiclassical limit of 
the Maxwell equations |j] . Does a similar procedure work 
for our model ? The question is open. 

Our theory is neither relativistic nor non-relativistic, 
since it does not involve time at all; it is based on the 
Euclidean group - which is indeed a subgroup of both 
the Galilei and of the Poincare groups. Our Euclidean 
model arises in fact as a reduction by time translations of 
the zero-mass spinning orbits of both the Galilei and the 
Poincare groups [Tlj - which constitute the conventional 
descriptions of "classical light" p^Hof). 

An application of our semiclassical model is the deriva- 
tion of the modified laws of refraction and reflection at 
the interface of two homogeneous media with different re- 
fractive indices. As found by Onoda et al. [§], polarized 
light suffers, in fact, a transverse shift. This "Optical 
Hall effect" [3 is indeed an optical version of the spin- 
Hall effect [3. Their shift formula can be rederived [ll|, 
following Souriau 0] , who argues that the two "mechan- 
ical" states on both sides of the interface are related by a 
symplectic transformation, S, which is indeed the classi- 
cal counterpart of the quantum scattering matrix. This 
transformation commutes with the symmetries of the op- 
tical device; in our case, this is plainly the Euclidean 
group generated by translations of the separating plane 
and rotations around its normal direction, N. Tedious 
calculation provides us with the "classical scattering ma- 



trix", s H3- 

Firstly, the conservation of planar linear momentum 
extends Snel's laws to spinning light, namely 



n- m suit 



= n out smtf out , 



(25) 



for refraction, and 9- m = ir — out , for reflection, respec- 
tively, where ni n resp.d n out denote the refractive indices 
on both sides of the interface. 

Next, the conservation of the (planar) angular momen- 
tum implies that light is shifted transversally 0,0], viz. 



[s out COS 6» ou t - S in COS in ] N X U in 
q OU t - qin = 1 — „ 1 itvt.... I ' ( 26 ) 



k m n sm ( 



IN x Ui , 



where Si n = s ou t for refraction. 

Notice that the shift depends in general on the wave- 
length, A = 2ttH/k, white light is split, in general, into 
colors, shifted by different amounts. For n out = — n- m , 
however, Snel's laws entail that the shift i|26|) vanishes: 
white light is not decomposed and is indeed refracted 
following the classic Snel law, as if it had no spin ! 

This case is not of pure academic interest, owing to the 
existence of left-handed media (with a negative refractive 
index) In the ideal case, one can have n = — 1, and 
a simple slab with parallel sides provides us with a 
"perfect lens" with no chromatic aberration. 

We note that the shift 1)26(1 vanishes also for a reflec- 
tion, since then s, n = — s out . This does not contradict the 
results in Imbert |l5l |. which are indeed of higher-order. 
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